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Abstract 

A quantum field theoretic treatment of inclusive deep-inelastic diffractive scattering is 
given. The process can be described in the general framework of non-forward scattering 
processes using the light-cone expansion in the generalized Bjorken region. Target mass 
and finite t corrections of the diffractive hadronic tensor are derived at the level of the 
twist-2 contributions both for the unpolarized and the polarized case. They modify the 
expressions contributing in the limit t, M 2 — > for larger values of (3 or/and t in the region 
of low Q 2 . The different diffractive structure functions are expressed through integrals over 
the relative momentum of non-perturbative t-dependent 2-particle distribution functions. 
In the limit t, M 2 — > these distribution functions are the diffractive parton distribution. 
Relations between the different diffractive structure functions are derived. 
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1 Introduction 



Deep inelastic diffractive lepton-nucleon scattering was observed at the electron-proton collider 
HERA some years ago [1]. This process is measured in detail by now [2] and the structure 
function F®(x,Q 2 ) was extracted. 1 The experimental measurements clearly showed that the 
scaling violations of the deep-inelastic and the diffractive structure functions in the deep-inelastic 
regime, after an appropriate change of kinematic variables, are the same. Furthermore the ratio 
of the two quantities, did not vary strongly, cf. [4]. While the former property is clearly of 
perturbative nature, the latter is of non-perturbative origin. 

The process of deep-inelastic diffractive scattering was first described phenomenologically 
[5]. Diffractive events are characterized by a rapidity gap between the diffractive nucleon and 
the remaining part of the produced hadrons, which is sufficiently large. Actually it is this 
experimental signature along with factorization for the twist -2 contributions [6] for this process, 
which allows to give a consistent field theoretic description. Due to this phenomenological 
considerations containing reference to specific pomeron models can be thoroughly avoided. In 
the limit of vanishing target masses the scattering cross sections and relations between the 
diffractive structure functions were derived in Refs. [7-9] for unpolarized and polarized diffractive 
scattering. In [7] we showed, that the scaling violations in the deep-inelastic and deeply inelastic 
diffractive case have to be the same due to the fact that the scaling violations are actually those of 
the operators which remain taking the respective matrix elements confirming the experimental 
observation. The set of structure functions which emerge in both scattering cross sections is 
actually larger than measured in current experiments. At low scales of Q 2 target mass effects 
become relevant similar to the case of deep-inelastic scattering [10-14], see also [15]. 

In the present paper we extend the picture developed in the massless case [7, 8] to the case 
of finite target masses and finite values of t on the level of twist-2 operators in the light cone 
expansion [16]. While in absence of mass effects the two-particle problem could effectively be 
reduced to a single particle description for the case t — * 0, this is no longer the case for finite values 
of t and/or target masses. Here two-particle effects become relevant, which do not allow for a 
direct partonic description. The variables t = (pi —p/) 2 and M 2 = p 2 = p 2 are closely connected 
and the simplification emerges if these scales vanish. Yet one may still follow the field theoretic 
picture developed in [7,8] in the general case M 2 ,t ^ and derive expressions for the diffractive 
structure functions including relations between them. At low scales Q 2 and large values of (3 
target mass corrections have to be considered in the experimental analysis. This generally applies 
also to finite values of t, unless the scale Q 2 is large enough. The diffractive structure functions are 
found as integrals over two-particle correlation functions f(z+, Z-\ t) between the incoming and 
outgoing nucleon. Here, z± denote the corresponding collinear light-cone momentum fractions 
and t is the relative momentum transfer squared between the incoming and outgoing proton 
momentum. We refer to the formalism of non-forward Compton scattering, cf. [17], and apply 
the general group theoretical algorithm of decomposing off-cone tensor operators into operators 
of definite geometric twist [18-20] to determine the contributions at twist 2. The analysis can 
be generalized to operators of higher twist. On the level of the various correlation functions 
relations can be established. These relations correspond to relations between structure functions, 
cf. [12,21-24]. In the case of deeply virtual Compton scattering relations of this type were found 
in [25-27] before. 

The present field-theoretic formalism for deep-inelastic diffractive scattering was also devel- 

lr The measurement of the longitudinal diffractive structure function F®(x, Q 2 ) has not yet been possible, but 
would be important. For the DIS structure function cf. [3]. Likewise the polarized diffractive structure functions 
gP 2 (x, Q 2 ) should be measured in the future to reveal the effects of nucleon polarization in this process. 
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oped in view of possible future measurements of the respective operator-matrixelements using 
lattice techniques, as successfully applied in case of the moments of the deep-inelastic structure 
functions [28]. We regard it as a challenge in future investigations to verify the experimentally 
observed ratio of moments in the diffractive and deep-inelastic case with these technologies. This 
ratio still awaits a rigorous non-perturbative explanation. 

The paper is organized as follows. In section 2 we describe the basic formalism. In section 3 
the symmetric part of the Compton amplitude is dealt with, through which the diffractive 
structure functions for unpolarized nucleons are derived. The polarized structure functions are 
determined in section 4. In section 5 we derive relations between different structure functions 
and section 6 contains the conclusions. 



2 Basic Formalism 

In order to compute the twist-2 target mass and finite momentum transfer corrections in po- 
larized and unpolarized deep inelastic diffractive scattering we briefly recall the notations and 
conventions used in previous papers [7, 8] by two of the present authors. The process of deep- 
inelastic diffractive scattering belongs to the class of semi-inclusive processes and is described 
by the diagram in Figure 1. The differential scattering cross section for single-photon exchange 
is given by 

dW = \ ips<S> £ $ L " w " ' <"> 

spins 

Here s = {p\ + I) 2 is the cms energy squared of the process and M denotes the nucleon mass. 
The phase space dPS^ depends on five variables since one final state mass varies. We choose 
as basic variables 

x = 91 = _J^_ (22) 

Q2 + W 2_ M 2 2 q Pl ' y ' } 

the photon virtuality Q 2 = —q 2 , t = (j>2 — Pi) 2 the 4-momentum difference squared between 
incoming and outgoing nucleon, a variable describing the non-forwardness w.r.t. the incoming 
proton direction, 

Q 2 + Ml - 1 gp_ . . 

X * = Q* + W*-NP = -^ X (2 - 3) 
for M\ > t, and the angle $ between the lepton plane x i and the hadron plane p 1 x p 2 , 

(p, x l).(p, x p 9 ) 
cos$ = ^-i irr^ — , (2.4) 

\Pl X l\\Pl X P2\ 

where p± = p 2 ±Pi, W 2 = (pi + q) 2 denotes the hadronic mass squared and the diffractive mass 
squared is given by M\ = (q — pJ) 2 . 
The momenta p± obey 

p\ AM 2 

For later use we refer to the variables rj and j3 defined by 



(p+ p.) = 0, y -± = — - 1 . (2.5) 



QP- 

v = — 

QP-\ 



-x 



' 2 



x 



= ^- = -<l, (2-6) 
2 qp- x F 
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as well as to the transverse momentum variable 



k-=v + -— , (g7r_) = 0. (2.7) 
V 



The variables x, Xj>, /3 and 77 obey the inequalities 



< x < x w < 1, < x < (3 < 1, (2.8) 

2 2/5 1 x 

-00 < 1 - - < 1 - — = - < -1 < 77 < < . (2.9) 

X XT] 2 — X 




Figure 1: The virtual photon-hadron amplitude for diffractive ep scattering 

For the spin averaged cross section, the leptonic tensor is symmetric. Taking into account 
conservation of the electromagnetic current one obtains [7] 

w s w s w s 
w% = -&w; + pIpIj^ + p^pIj^ + [pIpI + pI.pI] • (2-10) 

Here and in the following we do not assume implicitly, that azimuthal integrals are performed as 
sometimes done in experiment. In the latter case the number of contributing structure function 
reduces. 

In the case of polarized nucleons we consider the initial state spin- vector S\ = S, S 2 = — M 2 , 
only and sum over the spin of the outgoing hadrons. One usually considers the longitudinal (||) 
and transverse (_L) spin projections choosing 

S|| = (0;0,0,0,M) (2.11) 
S ± = (0;cos7,sin 7 ,0)M . (2.12) 

Here 7 denotes the azimuthal angle. The antisymmetric part of the hadronic tensor was derived 
in [8] and is given by 

WS 



w a = 

[11/ 


1 \P\ l iP2v P\vP2[i\ £ PiP2qS 


+ 


i [P^ £ uS Pl q 


T 1 
— Plv £ ^Sp 1 q\ 


+ 


W a 

1 \P2iJ, £ vSpxq P2v £ fJ,Spiq\ 


+ 


■ r t 

1 \P\ii £ vSp2q 


T 1 
~~ Plu £ [iSp2q\ 


+ 


W a 

1 \P2f± £ vSp2q P2v £ liSpiq\ ^4 


+ 


■ \ T T 

1 [Pl[l £ up lP 2S ' 


T T 1 
" Plu £ [Lpxp 2 S\ 


+ 


W a 

■ \ T T T T 1 vv 7 

1 [P2[i £ up lP 2S P2u £ [ip lP 2S\ 


+ 




1 £ [ivqS 



M 4 

wz 

M 4 
Wi 

M 4 

wi 

M 2 



(2.13) 
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We will specify below which terms of this general structures contribute in case of deep-inelastic 
diffractive scattering. The kinematic factors above are constructed out of the four-vectors q, pi,p 2 
and S as well as g^ u and e VQV1V2V3 using 

T 1-P T QfiQis lot 1 a\ 

P» = P„-q^, 9^ = 9^--^-, I 2 - 14 ) 

£ /AVlV2V3 ~ £ /J-VlV2V3 ~ £qviV2V'i~^ 1 (2.15) 

e TT = e - e ^ - e — (2 16) 

c /iuviV2 c ni/V!V2 c-qvviv? ^2 c t^qviV2^ 2 ' l i,IU J 

At the level of the twist-2 contributions factorization holds for diffractive scattering [6]. 
A. Mueller's generalized optical theorem [29] allows to move the final state proton into an initial 
state anti-proton, where both particle momenta are separated by t and form a 'quasi two- 
particle' state. The correctness of this procedure within the light-cone expansion relies, first, on 
the rapidity gap between the outgoing proton and the remainder hadronic part with invariant 
mass Mx and, second, on the special property of matrix elements of the light-cone operators 
which contain no absorptive part. The structure functions for the diffractive process can thus 
be obtained by analyzing the absorptive part of the expectation value 

Tp, (x) = (pi, -p 2 ,S;t\ % v {x) \p u -p 2 ,S;t) , (2.17) 



with T^ v defined as 



T^x) =iRT 



'!)*(-! ) 



:2.i8i 



As shown in [17,30] the operator is represented in lowest order of the non-local light-cone 
expansion by 



f„Jx) « -e 2 



' 2m 2 {x 2 -ie) 2 
where 



q aX n I ± _± \ _ aX n / _ __ 
Ofu/\ (J a \ 2' 2 / t iV I 2' 2 



(2.19) 



Sfj, v \aX — QfiaQuX + QfiXQua ~ 9^u9aX • (2.20) 

x denotes a light-like vector related to x, 

x = x — n[(nx) — \/ (nx) 2 — n 2 x 2 } , (2-21) 

with n a normalized time-like vector, n 2 = 1, and the bi-local light-ray operators O a and 
read 

O a (kiX, K 2 x) = i(fl a {k\X, K 2 x) - Sl a (k 2 x, kix)) (2.22) 

0\ (kiX, k 2 x) = £l 5 a (k,iX, k 2 x) + tt 5 a (k 2 x, Kix) , (2.23) 

with 

Q a (kix, k 2 x) = RT (: ip (kix) 7 q U (kiX, k 2 x) ip (k 2 x) : S) (2.24) 

Vt b a (k-iX, k 2 x) = RT (: ip (kix) 7 5 7 q U (ki£, k 2 x) ip (k 2 x) : S) , (2.25) 

where Ki = —k 2 = 1/2, cf. [17]. As is well-known, these operators contain contributions of up 
to twist four [18]. The scattering amplitude is obtained by the Fourier transform of the operator 
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Tim/} ( x ) an d forming the matrix element (|2.17j) . Here we want to study its twist-2 contributions 
including the target mass and finite momentum transfer corrections. This is obtained by har- 
monic extension of the twist-2 light-cone operators Q 1 ™ 2 (kiX, k 2 x) and fig™ 2 (kiX, k 2 x) to twist-2 
operators Q 1 ™ 2 {k\x, k 2 x) and fi^ 2 (kiX, k 2 x) defined off-cone. This procedure has been per- 
formed for QCD vector operators using group theoretical methods already in Ref. [31] and used 
for the target mass and t-corrections of virtual Compton scattering in Ref. [27], cf. also [32]. A 
general procedure for determining the complete (infinite) twist decomposition of off-cone vector 
operators was introduced in Ref. [19] and recently extended to arbitrary tensor operators in 
Ref. [20] . It should be emphasized that this method works at operator level before taking matrix 
elements. As a result the twist-2 Compton operator off the light-cone reads [27, 31] 

f^ 2 (q) = J d'xe^f^ix) (2.26) 

/djC g ^ X\ A 2 A 2 

2i7r 2 fx 2 ie) 2 v ^ ^ a ^ KX '' ~ KX ) ~ ^5a (kx, — KX)}, 



with 



f f d 

fi ($L («*> ~ KX ) = d cJ drj -^y 4 n( 5 ) M («) {^( 2 + xd ) ~ |i«r«V} (3 + xd)H 2 (u, ktx) , 

(2.27) 

cf. Eqs. (I2.22l2.23j) . and 

H v {u,Kx) = v 7 ^ (Ky/(ux) 2 - u 2 x 2 ^ l2 ~ V J v ^i 2 ^{ux) 2 - u 2 x 2 ^j e «M/ 2 , (2.28) 

and k = 1/2. In these relations a (formal) Fourier transform of the non-local operators is 
introduced and used off-cone later, i.e., 

{1 {5) ^kx,-kx)\x_ x = j c/ 4 ufi (5) »e iK ^ . (2.29) 

Let us emphasize first that (|2.27j) is an off-cone operator equation which determines the twist- 
2 part of fi( 5 ) M (/tx, — kx) by acting on Q tK i xu ) with a corresponding twist-2 projection operator 

Pa^ix.dx) whose result is presented above. Second, it has been shown in [31] that, due to this 
projection and the structure of the Fourier kernel, from the operator Q^^u) in (|2.27|) only its 
twist-2 part remains. 

To calculate the twist-2 part of the Compton amplitude T^ u (q) we have to parameterize the 
non-perturbative expectation values of the twist-2 operators 0^^(kx, —kx). Here, following 
the general approach of [31], this is done by requiring 

(Pi, -Vi\(? i (fi M ( M ) ^ bi> -Vi) 

= Z)*>(P±) / DZ6W(u-p-Z- -p + z+) f a (z-,z+,t) , (2.30) 

a J 

(pi, -p 2 , S\e 2 (fi 5(U («) + fi 5(U (— u)) \p u -p2, S) 

= Y, 1C *M> S ) / DZS^(u-p.z.-p + z + )f 5a (z.,z + ,t), (2.31) 

a J 

where, in fact, under the w-integration the generalized distribution amplitudes /( 5 ) a (z_, z + ; t) 
are reduced to their twist-2 contributions. These generalized distribution amplitudes depend 
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explicitly on t, and additionally there is a t— and M 2 — dependence which finally results through 
the distribution S^(u — p^z_ — p + z + ) from the Fourier transform in (|2.26|) . together with (|2.27l 
12.301 l2~3*T]) . cf. Section 3 and 4, as well as from the kinematic pre-factors /C" 5 ^(p±, S). 2 Here, 
for the symmetric part of the hadronic tensor we choose K, a ^(p±) as kinematic pre-factors 

K 1 " = pt, /C 2 " = ttI =K - V — . (2.32) 

V 

For the antisymmetric case the kinematic factors 

K\" = S», JC^=ffL^-, ^ = *-^T ( 2 -33) 



contribute. The normalization by M 2 in 1)2.33)1 is arbitrary and has to be arranged with 
the definition of the corresponding non-perturbative distribution functions. Furthermore, 
f(5)a(z-, z+',t) = f(5)a{^]t) denote the respective 2-particle amplitudes characterized by the 
two fractions z± of momenta p± and a relative kinematic separation in t. The momentum 
fractions of the incoming (z\) and outgoing (z 2 ) nucleon are formally unified into a 2-vector, 

Z = (z+, zJ) = ((z 2 + Zl )/2, (z 2 - Zl )/2) (2.34) 

with the measure 

DZ = 2 dz+dz- 9(1 - z + + z_) 9(1 + z + - z_) 9(1 - z + - z_) 9(1 + z + + z_) . (2.35) 

Due to the 5-distribution in (J2.301 12.31|) and the structure of the Compton amplitude, which is 
determined by ()2.26fJ2~2*%|) . it appears convenient also to unify the momenta as 

P = (p+,p-) = (P2 +PuP2 ~Pi) (2-36) 
and to abbreviate expressions according to 

if = «p*z = k^+^_) ) nj = y T ir = H/ , - , (2.37) 

which replaces km after performing the w-integration. Note that, due to the square roots in 
H.2(ku,x), the Fourier transform can only be performed if II 2 > since ^(qU) 2 - q 2 U 2 has to 
be real. In fact, this is equivalent to the requirement 

which restricts the allowed values of z± as a support condition to the generalized parton distri- 
butions /(5) (Z, t). It is convenient to introduce a common scale & for the integration variables 
z + = Z- = $(1 — C/v)> i- e -j to change them non-linearly into $ and £, 

0=*- + ^, C = ^- (2-39) 

Clearly, the \t\ values emerging in practice are not expected to deviate from this condition due to 
the presence of a non-perturbative damping factor exp(— b\t\) in the distribution functions, with 

2 In the following the explicit t-dependence of the distribution functions is understood also when we drop this 
variables in the respective expressions for brevity. 



7 



b ~ 4... 8 GeV -2 . Connected to the support condition the positivity of the variable II 2 (J2.47|) is 
guaranteed which is essential for performing Fourier transforms. The measure DZ* in these new 
variables reads 

DZ = 2 d$d(6(l - ■& + (1 + I/77) i?C)^(l + 1? — (1 + #C) 

6>(1 — — (1 — I/77) i?C)6>(l + 1? + (1 — I/77) ^C) • (2-40) 

The parameterization of the matrix elements (|2.3U|) and (|2.31|) is given by 

( Pl ,- P2 \e 2 i(n,(u)-Q,(u))\ Pl ,-p 2 ) = Y< Ka M) I DZ5(u-#V) f a (#,{;t)(2Al) 

( P i,-P2,s\e 2 (n 5fl (u) + n 5fl (u))\ Pl ,- P2 ,s) = ]T/q M (p ± ,s) fDZ6(u-#p) / 5o (AC;*)- 

a J 

(2.42) 

In the following we make the ^-dependence of the kinematics explicit : 

IP = K&V{rr,Q, (2.43) 

p»(j},0 = j£ + t£c=j£(i-CA?)+#C, (2-44) 

qU = K$(qV), (2.45) 

(qP) = (q P .)=q 2 /2(3<0, (2.46) 

n 2 = (k$) 2 P 2 , (2.47) 

p 2 = t(i - (/ry) 2 + (4M 2 - t) C 2 > 0, (2.48) 

(IT T ) 2 = (k$) 2 (P t ) 2 = (m?) 2 Wf-fP 2 (2.49) 

V 2 43 2 P 2 /Q 2 

P ' M > . (2.50) 



(pT)2 1 + 4/52 P 2/g 2 

"P 2 depends on both t and the integration variable ( which is limited by the Heaviside functions 
included in the measure ()2.4()j) . 

Recall, that V 2 > is required which for t = is trivially fulfilled. Concerning the region of 
small values of t < one observes that ( = is to be excluded. 

As will be shown in Sections 3 and 4, the absorptive part of the Compton amplitude solely 
emerges from a factor 1/[R(1) + ie] with the polynomial 

R(r) = (q + rU) 2 = (r 2 IT 2 + 2rqU + q 2 ) = U 2 (t - | + )(r - £_) . (2.51) 

The roots are 



- _ -gg ± ^(guy - g 2 n 2 _ q 2 

n 2 q u ± v/( g n) 2 - g 2 n 2 ' 



leading to 

1 1/1 1 



^(r) + ie 2 v /(glT) 2 - g 2 IT 2 \r - | + + ie r - £_ - ie 
The imaginary part of the Compton amplitude results thus from 

1 7T 1 



'2.53) 



Im 



R(l)+ie~ 2 y/(qU) 2 - g 2 n 2 L 
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+ . (2.54) 



Only the term containing the variable £ + contributes to the cross section, whereas there is no 
contribution due to the term This is seen as follows. For negative values of 1/rj < $ < 0, 
(qU) > 0, holds and therefore, in terms of variables (), one obtains 

L = ; 1 ^-^>0 (2.55) 

«0 1 + + A(3 2 V 2 /Q 2 K&- 

«0 1 - ^/T+WPVQ 2 k(3$V 2 ^ °° l ' j 

taking into account the support condition for V((). Here ==> denotes the limit Q 2 — > oo. 
Therefore 5(1 — £_) does not contribute. In the case of positive < $ < — 1/r/, one finds 

2/3 1 1 Q 2 

e+ = "^iT7iTlP/f ^^" + °° (2 ' 57) 

2(3 1 . . 

£- = ^ ; ~ < , 2.58 

K0 1 - ^1 + 4/? 2 P 2 /Q 2 ^ ~ 



and neither 5(1 — £ + ) nor 5(1 — £_) contribute. 
Having this in mind, let us rewrite 

i =i : 2/3 1 =— gvp2 (259) 

0' « 1 + + A(3 2 V 2 /Q 2 2k (3 1-^1 + 4^2^2 

from which it is evident that in case of diffractive scattering —2(3 takes the role of Bjorken 
variable Xsj in deep inelastic scattering as V 2 takes the role of AM 2 . Furthermore, rewriting 

5(i-£ + ) = hW-0 (2-60) 

it is obvious that after integration over d that variable is replaced by £ everywhere. For interme- 
diate momenta £ substitutes the scaling variable (3 and plays the role of a generalized Nachtmann 
variable. Let us note that this variable, however, implicitly depends on £, (3 and 77 through V(Q. 
Due to Eqs. (|2.52|) . the following equalities 

(H 2 + gn) 5(1 - £ + ) = y/(qU) 2 - q 2 U 2 5(1 - £ + ) = - (q 2 + gH) 5(1 - £ + ) (2.61) 

hold and can be used to simplify some expressions resulting for the absorptive part of the 
Compton amplitude. 



3 The Symmetric Part of the Amplitude 



In the case of unpolarized scattering the symmetric part of the Fourier transform of the operator 
x a 0^ 2 (kx, kx ) / (x 2 - ie) determines the contribution to the hadronic tensor of diffractive 
scattering. It reads (cf. [27], Eq. (6.1)), 

T^ } (q) = I — — — 2 x a ie 2 (^ 2 (KX,-Kx) -^ 2 (-kx,kx) 



\x L 



ie 



2ie 2 I dr (1 - r + r In r. 
o 

2 (a 2 ) 2 



dSi 



KT J \ KT 



[(q + u) 2 + ie\ 3 



9l(" T ? ) + 



u 2 q 2 



[(q + u) 2 + ie} 2 ^ v 



(3.1; 
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We have to form matrix elements (pi — P2\T^y(q)\pi, —P2) for this expression using the non- 
perturbative matrix elements (j2.41|) of the operator i(Q(u) — Q(—u)). This representation of the 
operator matrix element, based on the generalized parton distribution functions / a (Z), has to 
be inserted into the Fourier transform of the Compton operator. 

After a straightforward, but tedious calculation one obtains for the part which contains the 
factor l/(R(l) + is), see [27], Eq. (6.4), 



tw2 



w(<?) - y 



DZ< 



\ qU 



gll 

n 2 



n T n T 



(n 



{u T y 



qlC a ^ 



R(l)+ie 



(3.2) 



(n T ) z qU (ir 

In principle there occur also terms oc l/(i?(0) + ie) which, due to the overall factor q 2 , do not 
contribute to the imaginary part and are therefore omitted here. Obviously, the whole expres- 
sion obeys transversality proving gauge invariance of the complete amplitude. The functions 
F a (Z;g,IT) are 

nr 2 ! 2 



(gn) 2 - g 2 n 2 



(3.3) 



F 2 a (Z) = f a (Z) 



3n 2 (gn + n 2 



3 [II 



212 



(?n) 2 - 9 2 n 2 J„ t 2 ' \ tJ ( q ny-- q 2 iP 



Jo r Jo r l ^ Tr l J 



F 3 a (Z) 



3F a _ 2 / a (Z), 

2(gn) 



(3.4) 



(gn) 2 /z 

(gll) 2 -g 2 n 2 2 Vr 



n 2 ( g + n) 



(gll) 2 - g 2 n 2 J r 2 

( ? n) 2 - ? w[ (?n + ?2)r(Z) + (?n) r$ r (7 



r 



(gn) 2 - g 2 n 2 
n 2 



F 4 a (Z) 



F«(Z) 



3 ^3 + 
r 2 \r. 



n 2 



(gll) 2 - q 2 U 2 



(gll + g 2 ) f*(Z) + (gll) ^ ^ P(~ 



(3.5) 
,(3.6) 
(3.7) 



The imaginary part of the Compton amplitude is given by 



lmT { % 2 } (g) 



— 71 

x 



dd 



l/r, 



[(qV) 2 -q 2 V 2 \ x l'< 



s{& - 



qJCa 
qV 

qfca _ VlC a 
qV V 2 



(pry 



qV (P T ) 



(7?T)2 



(3.8) 
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We change now to the variables () and express the above integrals using the following new 
distribution amplitudes : 



$i 2 >(tf,C) 



.1/7/ 



/ dy 1 f a {y 1 ,Q = # / — - C), (3.9) 



1/7/ 



(3.11) 



In the following we consider the imaginary part of the Compton amplitude which forms the 
hadronic tensor for deep-inelastic diffractive scattering. As already stated above, the imaginary 
part follows from Im {l/(R + ie)}. The d integral is easily carried out 



lmT { % 2 } (g) 



-7T 



2k 2 7 [(gP) 2 -g 2 P 2 ]V2 



X 



qV 



\qP V 2 



(V 1 

' qV (pT)2 KsU'U 



(3.12) 



where the substitution ILj = was used. Here the functions F°-(Z;q, U) = 

Fi(£, C> <?; K CP)i i = 1, ■ ■ ■ ,5, are obtained from (|3.3|) - (|3.7|) using (|2.61|) in a non-trivial way : 



F?{t,Q = $a(£,C) + 
F?(£,() = $a(£,C) + 



[( g p)2_ g 2p2]l/2 

3kP 2 

[(gP)2 _ g 2p2]l/2 



K 2[-p2l2 

^0 + 7-^ K^ l^O (3.13) 



(gP)2 _ q 2 V 2 

Zk 2 [V 2 } 2 

(gP)2 _ g 2p2 



*«(£,0 (3-14) 



m 2 F (i t 

(qV) 2 -q 2 V 2 a2 \r^ 



+ 



2(gP) 



+ 



[(qV) 2 -q 2 V 2 } 1 / 2 y T 2 
K P 2 



F a4 (e,C) = 3F a3 (e,C) 



[( g p)2_ g 2p2]l/2 

2g 2 P 2 



K p 2 



[(<?P) 2 - q 2 V 2 ' 



F a5 ^C) 



2kV 2 



\{qV) 2 - q 2 V 2 \ x l 2 
Jo T Vr 



£$i 0) (£,c) + 



[(qV) 2 -q 2 
qV 

[(qV) 2 -q 2 V 2 \ x l 2 



qV 

[(qV) 2 - g 2 P 2 ]V2 



p2jl/2 V r ' ^ 



*i 1} (e,o 



(3.15) 



(3.16) 
(3.17) 
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In the limit t, M 2 — > 0, Eqs. (5,8) of Ref. [7], the hadronic tensor simplifies and contains 
only two structure functions which are linear combinations of the four occurring in the general 
case (I2.10J1 weighted by powers of (1 — xp). In this limit the functions F%\ k=3 do not vanish 
and contribute to the two structure functions remaining. On the other hand, F k \ k=3 do not 
contribute for forward scattering p 2 — > 0. 

Now, we can rewrite the imaginary part of the Compton amplitude, simply having recourse 
to [27], Eqs. (6.27) and (6.36), and replacing there x — > —2/3 : 



where 



ImT^ ( g ) 



lmT^ } (q) 



\mTlf u} (q) = \ m T^ } (q) 



lmT 2 ^ } (q) 



qV 



E 2?r / ^ 

a 

5>/«{( 



V 2 



qlCg VJCg 
qV V 2 



T 

9fiu 



-pT-pT 
o M " 
— T^|2 



(P 1 







X 






-h 




f i/ v 


+ 







rdiff 



(y,P,v,0 + 



Q 1 



A(3 2 V 2 



T T 
v H fia 



+ 



v T v T \ r a r 1 /^ dn 









v 2 J 



V£f(&fti?;C) 



-pT-pT 
T ' fi ' u 
9fa/ 



(pry 



(3.18) 



(3.19) 



(3.20) 



where 



and 



2/3 



v/1 + A(3 2 V 2 /Q 2 



2\3 



-2/3$i 0) (^Q 
+ 4(3 2 V 2 /Q 2 ' 

-2/5 2 P 2 9 / -2/3$i 1} (e;C) 



KT(e,/5,r/;C) 

Kf(£,/^;C) 



Q 2 dp \^l + A(3 2 V 2 IQ 2 



1 + 4(3 2 V 2 /Q< 



4(3 2 V 2 /Q 2 



(3.21) 
(3.22) 



:3.23) 



(3.24) 



Note that £ = f/;C) an d = V(r];Q. Differential operators as emerging in ()3.24|) are 
characteristic for target mass corrections and were applied earlier in Refs. [11,12], see also [23]. 
If the corresponding differential expressions are found the associated integral representations are 
derived straightforwardly. 
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We now convert the result referring to bases like used in (J2.1(J|) . By observing 



P^C = VJCa 1_ 

V 2 V 2 A(3 2 V 2 /Q 2 qV 



(3.25) 



one obtains 



± lm T { , u} (q) = E/^[- 9l () + ^ (£, o 



(3.26) 



T 3 a (e,o + ^T 4 a (e,o 



r T -P T + r T P T 



J?2 



(V T ) 2 \qV ' vs '^ V 2 
The distribution functions T?(£X) are related to the distributions used above by 
?i & = Wf? (£, /?, 77; C) - T 2 a (^, C) , 

' [2 w* ff - w*f + kT] (v, A m + o ^ & A ^ 



T 2 a (e,C) 



^o = W +w a d i ff )(e,/3,7 7; c)- 



Q 2 



2/3 2 P 2 



1A? ^ 



2 



(3.27) 
(3.28) 



3T 2 a (£,0~ Kf(£,/^;0 



T 4 a (e,0 = -2 



JW™(y,p,v;0+{3T^,()-V a d ?(t,p, m 



(3.29) 
(3.30) 

(3.31) 



The j3 and 77 dependence of Tj a (£, is understood. There are different possibilities to choose the 
tensor-basis of the hadronic tensor. One way is 



— lmT w (g) 



T 



^^/^^E[^^,0 4 



M 



VIC 



M 2 u 
1 



M 2 



M 2 



-jj, —V 



M 2 



' /V « rpal 



(P T ) 2 ^ L qV 



(3.32) 



Here the two kinematic invariants are K} p = p p _, K? p = ir p _. Introducing the structure functions 
Ui(j3,r]) one obtains 



2tt 



1 flU 



M 2 



„7T"„ + P"„7T"„ 7T T „7r T „ 

U 3 ((3,r})-> ~ v 



M 2 



M 2 



(3.33) 
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with 



d( 



(3.34) 



d( { Wf + 1ST Tl (e ' + 1*T T ' & 0] + ^ Tfa, } 



(3.35) 



U 3 (P,ri) 



M 5 



7?2 



M 2 



f/4(/3,77) 



M 2 



V 2 



Vtx 



(3.36) 



P 



M 2 1 



T 2 (£,C)J+^T 2 (£,C) 



(3.37) 



One may finally change to the basis in (|2.10J) . The unpolarized diffractive structure functions 
are then obtained by 



U 2 _ 2 ( 1 + ~ ) U 3 
U 2 + 2 ( 1 - ~ ) U 3 



1 + - u 4 

TjJ 

1-- u 4 



-U 2 + -U 3 
V 



1-^ 



(3.38) 
(3.39) 

(3.40) 
(3.41) 



with U = Ui({3,r}). The structure functions U, resp. W?, depend only on the measurable 
kinematic variables (3 and rj and form the hadronic tensor (J2.1(Jj) in the spin-averaged case for 
general azimuthal angle and can be measured directly. 

4 The Antisymmetric Part of the Amplitude 



The Fourier transform for the antisymmetric part of the amplitude again is taken over from 
Ref. [27] (cf. Eq. (5.2)): 



T [H (?) = \ E / DZ ^ 



-r 



[(gll) 2 - g 2 II 2 ] 



q a 1C a 5(3 F? 5]l (Z) + q a Il (3 (qJC a 5 )F? 5)2 (Z) 



with 



^i?(z) = (gn + n 2 ) / ^/ 5a F)+n 2 / ^ / ^ ha[ — 



1 dr „ /Z 



1 dr Z" 1 c?ri 



i?(l) + ie 



(4.1) 



o T Jo r i 



Z 



TTi 



(4.2) 
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/3gn(gn + n 2 



1 dr /Z 

2 / 5 a I 



/5a(z) ~~ V[( g n) 2 -g 2 n 

3 (gn) n 2 [ ld ^[ 1( ^±f (^_\ 

[(qU) 2 -q 2 U 2 } y r 2 7 r 2 MrrJ ' 



(4.3) 



^ ( ?(z) 



/ 5 „(z 
+ 



3gIT(g 2 + gIT) f 1 dr /Z 

2~ J5 a | — 



[(gn) 2 -g 2 n 2 
3 (gn) 2 



[(gn) 2 - g 2 n 2 



Let us remark that the symmetry (Z) — > — (Z) in i 7 ^ is a consequence of the symmetry in f 5a 
as has been shown in Ref. [27], App. C. Consequently, all terms oc 1/[R(0) + ie] drop out. 

We consider the imaginary part of this expression. Again, only the terms oc 1/[R(1) + ie] 
contribute, 



' m ^~[H (<?) = iv a/3 E/^/°0 2 

a J Jl/V 



(■&K) 3 [((fP) 2 - q 2 V 2 fl 2 



x 



(4.5) 



It is convenient to introduce the distribution amplitudes $5 a (*)(^, C) 
= 0/6.(0,0, 



1 dn 



'1 



.2 • • • 



1 dTi_i f 1 dn , 

2 J 5 a 



•1 ./0 



JO r l 



* yi 

1/v dy 



1 dTi-i 



V 
l/t? 



Ti • • • Ti_iTi 



r i-i Jti/i-n-n-!) Hi 



i) 



fba (Vi,0 



I/17 



i-l 



_ 2 y?-i y« 



l/j? 



fe/t/5o(j/t,C)] 



(4.6) 



C _1) (y,C), 

y 

where a successive change of variables y m — i?/(ti . . . T m _iT m ), m = n, . . . , 1, has been made and 
the support restriction of / 5a (i?, £) to I/77 < $ < has been observed. 

Then, performing the ^-integration with the help of the 5-function, one gets the following 
expression which simply may be obtained from [27], Eqs. (5.10) replacing there x — > —2/3: 



"mia(?)=| c ^/ dC (f 

Qa ^5/3 



qV 



[1 + A(3 2 V 2 /Q 2 fl 2 

(1 - kv 2 /q 2 ) $2(e, - /^7q 2 



(4.7) 



+ 



q a Vp{qK 



+3 



gP qV 

^V 2 /Q 2 



[l + 4/5 2 P 2 /g 2 ]V2 



*£ } (£,c) 



[1 + A(3 2 V 2 /Q 2 ] 



(4- 
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-ft 



+ 



V 2 q a Vp (VJC a 5 ] 
Q 2 qV V 2 

2 - Af3 2 V 2 /Q 2 _ w 



[i + 4/3 2 p 2 /g2]i/2 



[i + 4/5 2 P 2 /g 2 ] 

This equation is written more conveniently as 3 



-7r e 



/(!' 



^C<|^^^i(/3,C) + ^ 2 (AC)] 



2Q 2 P 2 9a0{P ^ ) 



qV qV 



<fe(/?,C) 



Here the following equalities 



1-P£P 2 /Q 2 

4 *<<>«; 



a/? V £ [i + 4/3 2 P 2 /g 2 ]V2 



(4.9) 

(4.10) 

(4.11) 
y/l + 4/3 2 P 2 /Q 2 

-^TTWPW n °' <412) 

2/3 1 - [1 + 4/3 2 P 2 /Q 2 ] 1 / 2 _ 1 Af3 2 V 2 /Q 2 
£ [1 + 4/3 2 P 2 /Q 2 ] 3 / 2 ~ 2 [1 + 4/3 2 P 2 /Q 2 ] 3 / 2 ' 



'1 + 4/3 2 P 2 /Q 2 = -(1 + 4/5/0 



were used. For the functions g a i(P, C) one finds the following expressions: 



-P 



d_ 

dp 



d_ (2(3 sgfeQ 
P d/3 ^£ [l+4/? 2 P 2 /Q 2 ]V2 

2/5/e 



[i + 4/? 2 P 2 /g 2 ] 3 / 2 



[l + 4/3 2 P 2 /Q 2 ]V 2 



+ftV 2 /Q 2 



(2 - 4/3 2 P 2 /Q 2 ) 
[1 + Aj3 2 V 2 /Q 2 } 



(4.13) 



(4.14) 



[1 + 4/3 2 P 2 /Q 2 ] 3 / 2 
ftV 2 /Q 2 



1 + 2f3£V 2 /Q 2 
[l + A(3 2 V 2 /Q 2 ] 1 / 2 



+3 





[1 + /3 2 P 2 /Q 2 ] 
<9 ^2/3 



9/3 \ C [l + 4/3 2 P 2 /Q 2 ] 1 / 2 y 



-2/3/e 



[l + x 2 P 2 /Q 2 ] 3 / 2 L 



(i - ftv 2 iQ 2 ) C) - /^7Q 2 *£ } (£, 



(4.15) 



(4.16) 



5 The implicit ^-dependence has been suppressed for brevity. 
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SaO (AC) = P 



d^ 2 l 4/? [l + 4/32p2/g2]l/2 

4/3 2 



[1 + 4/3 2 P 2 /Q 2 ] 3 / 2 
+ 



2 - 4/? 2 P 2 /Q 2 ^ (2) 



[1 + A(3 2 V 2 /Q 2 ] 



(4.17) 



Obviously, the contributions to the various kinematic expressions are related to only one function. 
Namely, introducing 



i [i + 4/? 2 P 2 /g 2 ] 1 /2 



(4.18) 



we may rewrite the functions g a i as follows 



d f„d 



(4.19) 



9az(P,V,Q = 

9ao(P,V,Q = 

One finally has to perform the (^-integral to obtain the diffractive structure functions in the 

2 1 2 

polarized case. The distribution functions g a i{P,7],Q\ i=0 occur in products with C L_ an d 
result into six structure functions of the form for each invariant a 



G aj (P,v) 
K aj (P,v) 



d(g a j{P,v,0, 

d( C9aj(P,ViO, 

d( C 2 9aj(P,V,0i 



(4.20) 
(4.21) 
(4.22) 



leading to 



ImT^(g) 



-7r e 



I"' 



qp 

Q a p-p ( qK% 



qp^ \qp 

q a 7T— /3 ( q)Q 
qp_ 



Ga2(P,V)~ 



2 Q 2 



2 g 2 



qp~ 
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Inserting the three invariants /C a 5 explicitly one obtains 



+ 



QaP-0 

qp~ 



+ 



2Q 2 



qp^ 
qS 



qp- 



G 12 + 



P2S 
M 2 



G21 



+ 7^ [pLG 20 + P-7T-(G 30 + H 20 ) + 7r!#3o) 



QS rr fa \ , 



M 2 G 10 + M 2 ^— ^ 
rj 



G31 + G% 2 — H 22 



p 2 _H 20 + p_7T_(#30 + ^20) + TT^SO + M 2 i/ 1( 



,77 - 1 



# 10 



(4.23) 



with 



P- 



P-TT- 



-t/v, 



7T 



4M 2 - \jif 



77 - 1 _ 2/3 
77 x 



(4.24) 



This structure of the hadronic tensor specifies the general structure obtained for three charac- 
teristic 4-vectors P\,p 2 and S" in Eq. (J2.13|) for the present process. One may cast (|4.23|) into 
the form of (j2.13|) using Shouten identities [33]: 



A ~fj,£ v par 

QXp^VpOT 



X v £ fj,pcTT p& V '/ICTT -^-a^-uppr -^T^-upap 

+ g\r 



(4.25) 
(4.26) 



Note that the structure functions Jko, J = G, H, K are suppressed by fi 2 /Q 2 relative to the 
other contributions, with fj, 2 a hadronic scale. In the Bjorken limit, these structure functions do 
not contribute. 



5 Relations between Diffractive Structure Functions 

Relations between different structure functions as observed in deep-inelastic scattering [12,21-23, 
25] hold also for polarized and unpolarized diffractive scattering in the limit M 2 ,t — > 0, cf. [7,8]. 
The situation is more involved for the case studied in the present paper, since the structure 
functions emerge as a (^—integral of sub-system structure functions, which accounts for the two- 
particle nature of the wave-function. Thus the corresponding relations can be established for 
the un-integrated (^-dependent functions only. 

5.1 Unpolarized Case 

As before in the case of deep-inelastic scattering [11, 12] and non-forward scattering [27] we seek 
for a representation of the distribution functions W^ L , V^S ^ in terms of a single function. One 

may choose the distribution functions $a 2 ^(£;C) (|3.10|) and express the functions W^2 L , ± \ 
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through the following differential relations : 



W^(£(3) Br,- C) ~ - 2 ~ (~ W ) (5 I) 

= ?U£* f«™o\ ( , 2) 



Note that the longitudinal distribution function W^(^(f3), (3, rj; £) vanishes in the limit 

p 2 /g 2 -> o. 

A generalized Callan-Gross relation, which holds for diffractive scattering in the limit 
M 2 , t — > 0, [7], is broken as in the case of deep-inelastic scattering [11]. The distribution 
functions WTjTj) are related by 

wffm,P,mO + w£m,P,mO = (1 + ^ m ^ wff (flfl, A m 



n 2 

p2 



wH?(t(j3),P,v,Q • (5.3) 



The distribution functions V^S ^ can be expressed by the distribution functions Wf^ff directly 



i (0.1) 



1 - ; ) Wai ff (C(/3),/9,^C) 

TTWvyo 2 



4(3 2 V 2 /Q 2 f l dp 



[1 + 4^2p2/g 2 ]3/2 



/ %W^(^p),(3p, V ;0. (5.5) 

J 13 P 



Three of the above distribution functions are independent. To obtain the four diffractive struc- 
ture functions on the level of observables iy/|f =1 in the unpolarized case (j2.10j) the (^-integral 
has to be performed. As shown in section |3] the respective linear combinations are weighted 
by different (^-dependent functions, that in general no relations exist on the level of structure 
functions. 



5.2 Polarized Case 

Below the (^—integral the distribution functions g a i\ 2 =o are all functions of T a (P,r],C) (EH). 
Therefore one may express the functions g a o and g a2 in terms of g a \ as the central function. This 
is achieved applying the relations 

P-OsW'V'Q = - T-^f^C), (5.6) 

op Jp y 

Fa(P,v,C) = - f- f d 4^\y\vX). (5.7) 
jp y J y y 

We consider twist-2 operators only and find the Wandzura-Wilczek relation 

9% 2 (P,vX) = -g t I 2 ((3,vX)+ I' -g% 2 (y,vX) (5.8) 

h y 
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between g^ 2 and g^ 2 ■ All target mass and t-corrections are absorbed into the structure func- 
tions. Note that this relation holds for all invariants a independently. The validity of the 
Wandzura-Wilczek relation also in case of diffractive scattering at general hadronic scales M 2 , t 
is a further example in a long list of cases: covariant parton model [23,24,34], target- and quark- 
mass corrections [12,14] gluon-induced heavy flavor production [35], non-forward scattering [25], 
diffractive scattering in the limit M 2 , t — > [8]. 
The distribution function g^ 2 obeys the relation 

9%*(P,v,0 = -m&iP^Q - ji ^ 2 v 2 /q 2 Y/ 2 jf iffi 2 ^' 

[1+4/3^/^3/2^ y J y yl 9ai (V,V,Q- 

It contributes to parts of the hadronic tensor which are suppressed by an overall factor fi 2 /Q 2 , 
with /i a hadronic mass scale. In the related case of non-forward scattering the emergence of 
this function has been observed in Ref. [25]. 

For g a 2 the ^-integral can be performed maintaining the structure of relation ()5.8j) . Therefore 
the Wandzura-Wilczek relation also holds for the diffractive structure functions 

G&(P,V) = ~GT(M+ C—Gfffav), (5-10) 

j/3 y 
j/3 y 

Due to the ^-dependence of g a o through £ and V 2 one needs to know g a i(P, Vi to calculate the 
structure functions J a> o([3,r]). As the integral (j4.20J) cannot be inverted, relations between the 
structure functions J a> o and the other polarized structure functions cannot be established unless 
referring to the un- integrated distribution functions g a i({3,T),C)- 



6 Conclusions 

Deep-inelastic diffractive scattering, like other hard scattering processes off nucleons, requires 
target mass corrections in the region of lower scales of Q 2 . In fact the nucleon mass M is not the 
only hadronic scale relevant to this process in which both the incoming and outgoing nucleon 
play a role. The modulus of the invariant t = (p2 — Pi) 2 on average is of the same size as M 2 . 4 
It appears to be natural to seek for a formalism dealing with both the target mass and t-effects 
simultaneously. Diffractive factorization and the use of A. Mueller's generalized optical theorem 
made it possible to reformulate diffractive scattering in terms of deep-inelastic scattering off 
an effective two-nucleon state accounting for t. We applied the formalism of the (non-local) 
light-cone expansion in the case of non-forward scattering and adapted it to the kinematics 
present in deep-inelastic diffractive scattering. The formalism accounting for target mass and 
t-corrections was then implemented generalizing the picture derived in the limit M 2 ,t — > in 
earlier work both for the unpolarized and polarized case. The present formalism can in principle 
be extended to the case of higher twist concerning the corresponding operator expressions f^. 
However, besides factorizable contributions also non-factorizable terms have to be described. 

In case of related semi-exclusive processes in which more than one final-state hadron is well separated in 
rapidity from the inclusively produced hadrons other invariants more would emerge. 
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While in the limit M 2 ,t — > all diffractive structure functions can be expressed in terms of 
diffractive parton densities, this is no longer the case at lower scales of Q 2 in regions where M 2 
and t-corrections become relevant. Here we mean the kinematic t contributions which have to 
be distinguished from the dynamic contributions emerging in the non-perturbative distribution 
functions. The reason for this is that the 2-particle kinematics of the incoming and outgoing 
nucleon deviates from quasi-collinearity which was instrumental to unify the boost variables 
of both particles to a single one. The condition for the absorptive part is now related to a 
5-distribution with a more complicated argument leaving the integral over the relative trans- 
verse motion of the two nucleons. This integral cannot be inverted in general referring to the 
accessible kinematic variables. The corresponding distribution functions 3> a (5)(£> T), £; t) have 
to be determined by non-perturbative methods in the future. 

Although the partonic picture does not hold for the diffractive structure functions one still 
finds it below the ^-integral. At the level of twist-2 the structure functions can be built from the 
corresponding operator expectation values as in the case of deep-inelastic scattering since the 
specifics of diffractive scattering is moved into the corresponding two-particle wave-functions. 
Consequently, the logarithmic scaling violations, which can be completely associated with that 
of the operators, cf. [7,17], are found to be the same as in diffractive scattering changing the 
variables x in the deep-inelastic case to (3 in the diffractive case. 

The presence of target mass and t-effects enlarges the number of structure functions deter- 
mining the hadronic tensor. In the unpolarized case four structure functions contribute, which 
cannot be related to each other directly. The Callan-Gross relation does not hold. However, in 
the polarized case the Wandzura-Wilzcek relation remains unbroken and holds even separately 
for the contributions to the different invariants A^ a |^ =3 . In addition to the structure functions 
surviving in the limit M 2 ,t — > 0, Q 2 — > oo several new structure functions contribute. They are 
damped however oc 1/Q 2 w.r.t. the other contributions. The present formalism can be used in 
experimental analysis of deep-inelastic diffractive scattering data referring to suitable models for 
the un-integrated distribution functions depending on (, for which rigorous determination using 
methods of non-perturbative QCD do not yet exits. In this way the structures being derived in 
the present paper can be tested. 

Acknowledgments: This work was supported in part by DFG Sonderforschungsbereich Tran- 
sregio 9, Computergestiitzte Theoretische Physik. 



21 



References 



[1] M. Derrick et al. [ZEUS Collaboration], Phys. Lett. B 315 (1993) 481; 
T. Ahmed et al. [HI Collaboration], Nucl. Phys. B 429 (1994) 477. 



[2] M. Derrick et al. [ZEUS Collaboration], Z. Phys. C 70 (1996) 391 |arXiv:hep- ex/9602010 ; 
C. Adloff et al. [HI Collaboration], Z. Phys. C 76 (1997) 613 |arXiv:hep-ex/9708016| ; " 
S. Chekanov et al [ZEUS Collaboration], Nucl. Phys. B 713 (2005) 3 
|arXiv:hep-ex/050106"0| . 

[3] C. Adloff et al. [HI Collaboration], Phys. Lett. B 393 (1997) 452 |arXiv:hep-ex/9611017| . 

[4] H. Abramowicz and J. B. Dainton, J. Phys. G 22 (1996) 911. 

[5] G. Ingelman and P. E. Schlein, Phys. Lett. B 152 (1985) 256; 

J. Bartels and G. Ingelman, Phys. Lett. B 235 (1990) 175; 

A. Donnachie and P. V. Landshoff, Phys. Lett. B 518 (2001) 63 larXiv:hep-ph/0105088| ; 
Phy s. Lett. B 191 (1987) 30 9 [Erratum-ibid. B 198 (1987) 590]; Phys. Lett. B 437 (1998) 
408 |arXiv:hep-ph/9806344| ; Nucl. Phys. B 244 (1984) 322; refe rences quoted therein; 
W. Buchmiiller and A. Hebecker, Phys. Lett. B 355 (1995) 573 |j arXiv:hep-ph/9504374| ; 
M. Wiisthoff and A. D. Martin, J. Phys. G 25 (1999) R309 |arXiv:hep-ph/9909362| ; " 
A. Hebecker, Phys. Rep. 331 (2000) 1 

J. Bartels, J. R. Ellis H. Kowalski and M. Wiisthoff, Eur. Phys. J. C 7 (1999) 443 
arXiv:hep-ph/9803497| ; 



J. Bartels and H. Kowalski, Eu r. Phys. J. C 1 9 (2001) 693 |arXiv:hep-ph700 10345 1 ; 
H. Navelet and R. Peschanski, arXiv:hep-ph/0 105030 and references quoted therein; 

G. Ingelman, |arXiv:hep-ph /9912534| and references therein; 

M. Diehl, T. Feldmann, R. Jakob, and P. Kroll, Phys. Lett. B460 (1999) 204; Eur. Phys. 
J. C8 (1999) 409; 

P. Kroll, Nucl. Phys. A666 (2000) 3. 

[6] A. Berera and D. E. Soper, Phys Rev. D 53 (1996) 6162 |arXiv:hep-ph/9509239| ; 
A. Berera, |arXiv:hep-ph/9 606448 

J. C. Collins, Phys. Rev. D 57 (1998) 3051 [Erratum-ibid. D 61 (2000) 019902] 
|arXiv:hep-ph/9709499| . 

J. Blumlein and D. Robaschik, Phys. Lett. B 517 (2001) 222 [arXiv:hep-ph/0106037| . 
J. Blumlein and D. Robaschik, Phys. Rev. D 65 (2002) 096002 |arXiv:hep-ph/0202077| . 
J. Blumlein and D. Robaschik, Nucl. Phys. A 711 (2002) 228 |arXiv:hep-ph/0207250| . 
O. Nachtmann, Nucl. Phys. B 63 (1973) 237. 

H. Georgi and H. D. Politzer, Phys. Rev. D 14 (1976) 1829. 



[7 
[8 
[9 
[10 

[11 
[12 

[13 



J. Blumlein and A. Tkabladze, Nucl. Phys. B 553 (1999) 427 |arXiv:hep-ph/9812478| . 

J. Blumlein and A . Tkabladze, Nucl. Phys. Proc. Suppl. 79 (1999) 541 
|arXiv:hep-ph/9905524| . 



[14] A. Piccione and G. Ridolfi, Nucl. Phys. B 513 (1998) 301 |arXiv:hep-ph/9707478| . 



22 



[15] V. Baluni and E. Eichten, Phys. Rev. D 14 (1976) 3045; Phys. Rev. Lett. 37 (1976) 1181; 
S. Wandzura, Nucl. Phys. B 122 (1977) 412; 
S. Matsuda and T. Uematsu, Nucl. Phys. B 168 (1980) 181; 

R. L. Jaffe and M. Soldate, Phys. Rev. D 26 (1982) 49; 

H. Kawamura and T. Uematsu, Phys. Lett. B 343 (1995) 346 [arXiv:hep-ph/9409455| . 

[16] S.A. Anikin and O.I. Zavialov, Ann. Phys. (N.Y.) 116 (1978) 135, 
M. Bordag and D. Robaschik, Nucl. Phys. B169 (1980) 445, 
M. Bordag, B. Geyer, J. Horejsi, and D. Robaschik, Z. Phys. C26 (1985) 275, 
O.I. Zavialov, Renormalized Quantum Field Theory, Dordrecht 1990, 

O.I. Zavialov, Theor. Math. Phys. 138 (2004) 370, [translated from: Teor. Mat. Fiz. 138 
(2004) 437]. 

[17] J. Blumlein, B. Ge yer and D. Robaschik, Nucl. Phys. B 560 (1999) 283 
|arXiv:hep-ph/99035"20| . 

[18] B. Geyer, M. Laz ar and D. Robaschik, Nucl. Phys. B 559 (1999) 339 
|arXiv:hep-th/9901090| ; 



B. Geyer and M. Lazar, Nucl. Phys. B 581 (2000) 341 |arXiv:hep-th700 03080 1 . 
[19] J. Eilers, B. Geyer and M. Lazar, Phys. Rev. D 69 (2004) 034015. 

[20] J. Eilers, The decomposition of local and nonlocal operators with repect to irreducible repre- 
sentations of the orthogonal group and some of its applications in Quantum Chromodynamics, 

PhD Thesis, Leipzig 2004; 

J. Eilers, Geometric twist decomposition off the light-cone for non-local QCD operators, to 

appear. 

[21] C. G. Callan and D. J. Gross, Phys. Rev. Lett. 22 (1969) 156. 
[22] S. Wandzura and F. Wilczek, Phys. Lett. B 72 (1977) 195. 

[23] J. Blumlein and N. Kochelev, Nucl. Phys. B 498 (1997) 285 |arXiv:hep-ph/9612318| . 

[24] J. Blumlein and N. Kochelev, Phys. Lett. B 381 (1996) 296 | |arXiv:hep-ph/9603397| . 

[25] J. Blumlein and D. Robaschik, Nucl. Phys. B 581 (2000) 449 |arXiv:hep-ph/0002071| . 

[26] B. Geyer and D. Robaschik, Phys. Rev. D 71 (2005) 054018 |arXiv:hep-ph/0407301| . 

[27] B. Geyer, D. Roba schik and J. Eilers, Nucl. Phys. B 704 (2005) 279 
|arXiv:hep-ph/04073"00| . 

[28] cf. S. Capitani, Acta Phys.Polon. B33 (2002) 3025; 
M. Guagnelli et al., Eur. Phys. J. A17 (2003) 365; 

M. Gockeler, R. Horsley, D. Pleiter, P. E. L. Rakow, A. Schafer, G. Schier- 
holz and W. Schroers [QCDSF Collaboration], Phys. Rev. Lett. 92 (2004) 042002 
|arXiv:hep-ph/03042"49| . 

P. Hagler, J. Negele, D. B. Renner, W. Schroers, T. Lippert and K. Schilling, [LHPC col- 
laboration], Phys. Rev. D 68 (2003) 034505 |arXiv:hep-lat/0304018| ; 
M. J. Savage, Talk given at 23rd International Symposium on Lattice Field Field: Lat- 
tice 2005, Trinity College, Dublin, Ireland, 25-30 Jul 2005, PoS LAT2005:020, 2006, 

|hep-lat/0509048l ; 

S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage, |arXiv:hep-lat/0602010 



23 



[29] A. H. Mueller, Phys. Rev. D 2 (1970) 2963; Phys. Rev. D 4 (1971) 150; 

P.D.P. Collins, An Introduction to Regge Theory and High Energy Physics, Cambridge 1977, 
pp. 331; 

A. A. Logunov, B. V. Medvedev, M. A. Mestvirishvili, V. P. Pavlov, M. K. Polivanov and 
A. D. Sukhanov, Theor. Math. Phys. 33 (1978) 935 [Teor. Mat. Fiz. 33 (1977) 149]. 

J. Bliimlein, J. Eilers, B. Geyer and D. Robaschik, Phys. Rev. D 65 (2002) 054029 
|arXiv:hep-ph/0108095 . 



[31] B. Geyer, M . Lazar and D. Robas chik, Nucl. Phys. B 618 (2001) 99 [Erratum-ibid. B 652 
(2003) 408] |arXiv:hep-ph/0108061| . 

[32] A. V. Belitsky and D. Miiller, Phys. Lett. B 507 (2001) 173 |arXiv:hep-ph/0102224| . 

[33] M. J. G. Veltman, Nucl. Phys. B 319 (1989) 253; 

B. De Wit and J. Smith, Field Theory in Particle Physics, (North Holland, Amsterdam, 
1986), p. 414. 



[34] R. G. Roberts and G. G. Ross, Phys. Lett. B 373 (1996) 235 |arXiv:hep-ph7 9601235 



J. Bliimlein, V. Ravin dran and W. L. van Neerven, Phys. Rev. D 68 (2003) 114004 
[arXiv:hep-ph/0304292| . 



21 



